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REMARK ON THE THEORY OF LIFTING SURFACES* 
By Aldo Muggia 

SUMMARY 


First of all, a brief synopsis of the Weissinger method, as it applies 
to a rectangular wing, is set forth, in order to show how lifting surface 
theory is applied in this simple case and to show that his idealization 
of the vortex system is justifiable in this particular instance, Ry 
building on this framework and merely adding a few approximations and 
unrestrictive understandings, it is demonstrated how the same sort of 
vortex system can be devised, and can find sanction, for the treatment 
of the lift problem presented by any thin wing of arbitrary plan form. 

1. To begin with, let attention be directed to the aerodynamics of 
a lifting surface (which is a suitable idealization of an actual wing) 
having the simple physical property that it departs but slightly from the 
flat surface S, which is the projection of the lifting surface on the 
xy -plane. Furthermore, this surface is to be considered Immersed in an 
incompressible fluid of density p and to have a free-stream pressure and 
velocity denoted by p^ and V., respectively. This impinging stream 

is assumed to be directed along and have the same sense as the positive 
x-axis. In the right-handed Cartesian coordinate system employed here 
(see fig. 1 ) it will be assumed that the z-axis is in the direction of 
the vertically downward pointing vector, as illustrated. 

On the basis of the above-stated hypotheses, it is legitimate to 
assume that the perturbations to the free-stream uniform velocity V m 
that are produced by the presence of the lifting surface, will be small. 
Consequently, it follows that the local pressure p will be an harmonic 
function (a solution to Laplace's equation) of the position-coordinates x, 
y, z, and, in addition, the overpressure at such a point will be dependent 
upon the potential cp describing the behavior of the local incremental 
velocities, through means of the relation 


p - p 
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*"Sulla teoria delle superfici portanti." Atti della Accademia delle 
Scienze di Torino, vol. 87 , 1952-1953. Introduced by Carlo Ferrari, Active 
National Member, at the Session of 13 May 1953. 
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The pressure jump occasioned by passage from the underside to the 
topside of the surface will be a certain unknown function of the points 
of S; call it f(x,y). The value of this function must become zero along 
the edge of the surface S, and thus one may write that the overpressure 
is given by 


p(x,y,z) 


_ _1 p P f(x',y')z dx'dy' 
^°° tat J J g r 3 


where 


r = ^(x - x*) 2 + (y - y ') 2 + ^ 

Upon invoking the stipulation that one must have cp = 0 at an 
infinite distance upstream from the surface, it is seen that the sought 
potential will have the form 


<p(x,y,z) 



f (x' ,y ') dx'dy ' 

r^ 


Now let it be assumed, for convenience* sake, that the equation 
denoting the leading edge of the lifting surface is to be taken as 
x = x-^(y) and that for the trailing edge as x = X£( y), while, likewise, 

as is customary, the semispan is to be denoted by b. Then, upon 
carrying through a few obvious transformations it is possible to rewrite 
the expression giving the perturbation potential as 


cp(x,y,z) 



a,,, r*' 


f(x,y*)dx - 



p +0 ° dxi r x 2( y') 
Jx 2 ( y') r5 ^(y 1 ) 


f(x,y')dx 




Furthermore, let the following conventions be hypothesized: 

Take S' to be the wake region of the xy-planej ie., the region 
which lies downstream of the S region and lying in between the two 
straight boundaries y = -b. 
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Then make the definitions that on the wing region S 

, x 


m 


(x>y) = f f(x,y)dx 


while on the wake region S’ 


m 


( y ) f 

PVoo U -y 


X2(y) 


*i(y) 


f (x,y)dx 


Consequently, it may be seen that the perturbation potential is now 
expressible as 


cp(x,y,z) 


z r r m dx'dy* 
J S+S' 


which merely says that the sought value of cp is the potential which 
describes the behavior of a distribution of doublets of strength m 
per unit area in the S+S' region. 


This system of doublets may be replaced by a system of vortices, 
spread over the same S + S* region, by having recourse to the equiva- 
lency principle between doublets and vortices. The local strength of 
this equivalent vortex distribution, per unit distance, will have x- and 
y-components given by the following expressions: 

In the S region. 


7x 


_ ^m _ 


f: 
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dy pV m J X;L (y) $y 


dx 


and 


/y dx pV* 


while in the S* region. 


7 x = 
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The above-described system of vortices will induce, at any arbitrary 
general point of the S region, a certain velocity, the vertical component 
of which will be given by the integral 


v z 



£f I 

ay' y - 




Now this vertical component of the induced velocity must be of such a 
v z a z 

magnitude that — = holds true, in order that the velocity vector 
Voo dx 

representing the total flow at the point in question shall be tangent 
to the wing surface. Thus one now has obtained an integrodifferential 
equation to work with for the determination of the unknown function, 
f (x,y). 


2. The exact solution of this equation will, however, entail rather 
formidable difficulties, and for this reason it is best to fall back 
upon a much simpler approximate procedure for attaining the desired 
result. To this end, let the situation in regard to the simple rectan- 
gular wing be examined first of all (refs. 1 and 2). In this case one 
may let 

*1 (y) = -- and x 2 (y) = +| 

The angle of attack for any profile section situated at a spanwise dis- 
tance y from the center line of the wing may be denoted by a(y), where 
this symbol is meant to denote the true aerodynamic angle of attack of 
the profile in question, measured from the angle of attack for zero lift. 
Thus it follows that the governing integrodifferential equation may be 
written as 



Now make the approximation that the value of jx - x'|, that enters 
into the expression for r, is to be replaced by its average, which is 
simply l / 2 in this case. Making this substitution, and changing the 
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V / 




order of integration will result in the following simplified form for 
the integrodifferential equation: 


where 



4£_ . A y' 1 + 

dy' y - y' 2 



(y - y ') 2 


r(y') 


i r lf2 



ry^Sy')^' 


is the circulation function for the velocity distribution around the 
profile that is situated at the spanwise location denoted by y'. 

The above-derived approximate equation equates the vertical velocity 

component a(y)V 00 to the induced velocity produced at the point 

by action of a special system of vortices, which may be considered as 
derived from the actual distribution of vortices by concentration of all 
the bound vortices along the quarter-chord line. 

3. The line of reasoning followed in section 2 is only valid for 
the case of rectangular wings, but the result obtained is known to hold 
true even for other cases, as has already been pointed out by Weissinger. 
The wider generality of this result may be established upon the basis 
of the following considerations: 

Let it be assumed that the bound vortices are concentrated along a 
line denoted by x = xo(y), as illustrated in figure 2. The precise way 
in which this line is to be selected will be explained in full later on. 

It may be remarked here, however, that the acceptance of this represen- 
tation for the bound vortices is equivalent to replacement of the actual 
distribution of doublets in the S region by means of a special kind of 
doublet system, consisting of null doublets everywhere upstream of the 
x = xo(y) line and a distribution of doublets of strength (per unit 
area) described by means of the m(y) function throughout the region 

(and in the S' region as well), where the region is that area of 

the wing which lies downstream of the x = XQ(y) dividing line. The 

potential which describes this new redistribution of doublets is express- 
ible as 


<P*(x,y>z) = - 


4k 


II 


S-J+S 


m(y' idx'dy' 
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and thus the error in this approximate expression for the potential of 
the flow is given by the difference in the two potentials, i.e., this 
difference is 


<P(x,y,z) - «p*(x,y,z) = - 2 f f dx '? y f f(x,y')dx 

^jtpVoo J J g 

:IL 


4«pV„ 


dx'dy ' rW> 
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rx 2 ( y 1 ) 


px-2\ y ) 

/ f(x,y’)dx 

J x 1 (y*) 


po px ) p 

■r—rr- / dy' / f(x,y’)dx / 

^JtpVo oJ_ b , -'x- L (y) ^x 
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p» px 2 (y') px 2 i.y ' ) 

3-V / <*' / f(x,y')dx / flf 

‘tJtpVcou.'b Jx-^y’) ^ Xq (y 1 ) r 3 

_ C b , px 2 (y’) 


^x 2 (y 1 ) 


\/(x - x ’) 2 + (y - y')2 + z 2 


dx* 


x - x 0 (y') 


|/[x-x 0 (y)] 2 + (y -y ) 2 + z 2 Jx i^’) 


px 2 ( y ’ ) 

/ fCx'jy'Jdx' 

w v, M 


dy‘ 


Now let attention be focussed upon that region of space which is 
composed of all points which have quite small absolute values for the 
vertical coordinate |z| and which when projected upon the xy -plane 
fall within the S region; let this portion of space be labeled the 
2 control volume. For points within £, therefore, one may replace 
the cp function with the cp* approximation ^and thus it will be per- 
missible to substitute ^2— f 0r ^31 on -the surface of the wing^ pro- 

oz oz / 

vided the quantity standing within the curlicue brackets is of small 
enough size. 


Further, it is to be noted that for wings with sufficiently large 
r~"~" 2 2 2* 

spans, the value of y(x -x') + (y - y') +z does not vary to any 
marked degree as one ranges over the values of x' of interest, provided 

the distances ^(y - y ' )^ + z^ remain large enough, while on the other 

hand, if the distances ^(y - y ' ) 2 + z 2 are small, then the value 
that |x - x'| takes on in the Z control zone can be represented to 

good approximation by use of its average value ^ l (y) where the chord 
distribution function l (y) is defined as 

*(y) = x 2 (y) - x 1 (y) 
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Thus, in analogy to what was done in the case of the rectangular 
wing, it will he legitimate to make the approximation that 

^(x - x') 2 + (y - y') 2 + z 2 = /E " ^(y'H 2 + (y - y') 2 + z 2 

provided the point with coordinates (x,y,z) is so situated that it makes 
the relation 

|x - XQ(y)| = | l(y) 


hold true. 

If this is true, it follows, in consequence, that 

r x 2(y’) r n r x 2(y') 

/ f(x',y')(x - x')dx* = x - xq ( y • ) I / fCx^y'Jdx' 

J x 1 (y') L JJ x 1 (y l ) 

from which one obtains the desired definition for XQ(y') as 


r x 2 (y') r x 2(y') 

/ x'f (x',y')dx' / x'7 v (x',y' )dx' 

/ ,n ~ x i(y') ^ x i(y') 

*o ( y > * l - L " ■ ■■ • 


r x 2(y') r x 2W) 

/ ^ ( x ' ,y * )dx' / 7 y( x ' ,y 1 )^ x ’ 

^xiCy') ^ x i(y') 


' x 2(y') 


The interpretation of the relationship just deduced is as follows: 
The proper XQ(y') abscissa coordinate to choose at each profile section 

through the wing is the one which corresponds to the location of the 
barycenter of the circulation for that section. In other words, it is 
the barycenter of the moments of the vector forces fl taken about the 
points P(x’ > y'), where I is the unit vector in the direction of the 
x-axis and where P is the radius vector out to any arbitrary general 
point in the S region at which the circulation-function value is 
7 y ( x '>y'). 


Thus, to close approximation, one may select the Xq 
values according to the rule 


abscissae 


x o(y') = x i(y') + J- Hy') 
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while, when applying the boundary condition, it will be necessary, in 
addition, to make use of the potential values (or the induced velocity 
values) which appertain to the locations (x,y,0) for which it is true 
that 


x(y) = x Q (y) + | l(y) = XjJy) + | l(y) 

This result, which has been deduced by aid of the above-mentioned 
list of specific observations and series of approximations, may be arrived 
at by examination of the general equations applying to lifting surfaces. 
This result is important, for example, in those cases where one wishes 
to obtain the distribution of circulation (and thus of the lift) which 
exists out along the span of swept wings (refs . 3 and k ) . 


Translated by R. H. Cramer 
Cornell Aeronautical Laboratory, 


Inc . 
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Figure 2.- Location of the bound vortex line and areas of integration. 
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